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Summary 

Arrays  having  specified  be  am  widths  and  sidelobe  levels 
may  be  synthesized  from  a  simple -HSeed*^array  having  the  re-^N^  FIECTP  M 
quired  sidelobe  level,  but  a  much  broader  beamwidth  than  the  Ha  E“  H 

specification  value.  An  iterative  process  doubles  array  spaciiS^  DEC  1  0  198G9 
at  each  step.  This  narrows  the  main  beam  by  a  direct  scaling^T^N  jfl 

factoT  of  twain_sin^)  space.  The  echelon  lobes  introduced  by^V  Hj 

this  action  are  reduced  to  the  specification  sidelobe  level  values'1" 
by  an  interpolation  procedure  which  introduces  appropriate  ele¬ 
ment  weights  at  the  vacated  positions.  This  "stretch  and  fill" 
tactic  may  be  used  repeatedly  to  obtain  any  beamwidth  desired. 

The  computations  of  weights  for  these  arrays  of  any  size  are 
trivial  and  the  resulting  current  distributions  have  a  well- 
behaved  taper.  Aperture  size  requirements  and  efficiencies 
are  consistent  with  those  of  conventional  designs  for  a  given 
beamwidth  and  sidelobe  combination.  Supergaining  problems 
are  not  encountered  with  this  procedure. 

Interpolation  techniques  are  also  used  to  develop  a  unique 
aperture  format  and  excitation  matrix  combination.  This 
approach  to  beam  steering  provides  for  control  of  pointing  angle 
through  signal  level  control  as  opposed  to  the  conventional  use 
of  phase  shifters  for  steering. 
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SECTION  I 


BACKGROUND 

Since  the  work  presented  here  draws  heavily  upon  some  previous  interpolation  studies 
by  the  writer  [  1 J  ,  a  cursory  review  of  this  material  will  be  presented.  A  sampled-data 
(time-quantized)  time  series  is  considered  which  involves  a  sampling  period  T,  and  a  cor¬ 
responding  sampling  rate  R  (R  =  1/T).  The  power  density  spectrum  of  the  continuous  (but 
unavailable)  time  function  originally  sampled  is  assumed  to  be  band  limited  with  a  cutoff 
frequency  which  is  no  greater  that  the  foldover  frequency  (R/2)  of  the  sampling  process.  It 
is  desired  to  raise  the  sampling  rate  to  some  higher  value  by  some  efficient  means.  That 
is,  it  is  desired  to  estimate  the  sample  values  of  the  continuous  time  series  which  lie  be¬ 
tween  the  available  values  which  are  delivered  to  the  processor  at  T-second  intervals.  The 
estimation  errors  are  to  be  minimized  (in  some  sense)  for  a  given  data  processing  cost. 

The  procedure  employed  is  one  in  which  the  sampling  rate  is  doubled  at  each  of  a 
number  of  levels  using  processors  of  the  type  shown  in  Figure  1-1.  Time  samples  of  the 
continuous  function  s(t)  enter  a  tapped  delay  line,  and  the  weighted  sum  of  the  2N  samples 
is  used  as  the  estimate  of  the  intermediate  value  between  each  input  sample  pair.  The 
original  input  samples  are  recovered  at  the  a  tap,  delayed  one-half  the  input  sampling 
period,  and  merged  with  the  intermediate  estimates  from  the  summing  bus.  The  output 
sampled-data  stream  would  consequently  have  a  period  of  T/2  and  a  sampling  rate  2R. 

It  is  thus  argued  that  if  these  intermediate  value  estimates  are  the  best  that  can  be  obtained 
from  the  original  input  data,  then  these  estimates  should  be  used  if  further  sampling  rate 
multiplications  (additional  intermediate  estimates)  are  desired.  A  cascade  of  processors 
of  the  type  shown  in  Figure  1-1  are  indicated  with  sampling-rate  doubling  occurring  at 
each  level  or  stage. 


0  T  2T  t 

Figure  1-1.  Interpolation  Processor,  Sampled  Data  Version 


The  above  process  is  efficient  in  that  processing  results  obtained  from  previous  levels 
are  available  for  use  in  levels  following.  The  interpolator  design  at  each  level  can  be 
tailored  to  meet  the  needs  at  that  level.  As  the  sampling  rate  is  increased  through  the  cas¬ 
cade  the  interpolation  task  becomes  easier  since  the  available  input  samples  are  closer  to, 
and  hence  more  strongly  correlated  with,  the  intermediate  values  to  be  estimated.  Early 
in  the  cascade,  where  the  sampling  rate  is  low,  one  finds  the  the  most  complex  (longest) 
interpolation  processors.  The  simplest  (shortest)  interpolation  processors  are  found  at  the 
end  of  the  iterative  interpolation  chain  where  the  sampling  rate  is  highest.  Rate  multiplica¬ 
tion  factors  that  are  not  powers  of  two  may  be  obtained  by  easily-configured  multiplication/ 
resampling  procedures  (  1,  Section  5.4]  . 

The  performance  of  the  processor  of  Figure  1-1  may  be  evaluated  by  a  frequency- 
domain  technique  [  1,  Section  2.3]  which  leads,  concurrently,  to  a  method  for  determining 
the  processor  coefficients  of  Figure  1-1.  If  the  input  signal  s(t)  has  a  low-pass  spectrum 
which  extends  to  the  foldover  frequency  R/2,  the  spectrum  of  the  input  sampled-data  series 
will  contain  upper  and  lower  sideband  terms  of  total  bandwidth  R  located  about  suppressed 
(AM)  carriers  centered  at  positive  and  negative  multiples  of  the  input  sampling  rate  R.  If  a 
perfect  interpolation  could  be  achieved,  the  output  spectrum  would  be  identical  to  the  input 
spectrum  except  that  frequency  groups  at  odd  multiples  of  R  would  be  missing,  and  the 
groups  at  even  multiples  of  R  would  represent  a  doubling  of  (voltage)  level.  This  result 
follows  since  perfect  interpolation  would  be  equivalent  to  use  of  a  sampling  rate  of  2R  on 
the  original  s(t)  data. 

The  interpolation  portion  of  Figure  1-1  should  have  an  effective  (relative)  frequency 
response  (real,  even)  function  which  is  exactly  one  for  0  =£f  sr/2  and  is  minus  one 
from  R/2  ^  f  *  R/2.  A  good  design  approach  [  1,  Section  II]  employs  a  symmetric 
processor  of  total  length  2N.  That  is,  for  Figure  1-1 

a  -  a  (1-1) 

n  -n 

so  define 

b  2a  (1-2) 

which  yields  the  real  frequency  function  H  (f), 

N 

H(f)  =  £  bn  cos  (2n-l)  n  1  (1-3) 

n=l 
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Note  that 


H(f+2R)  =  H(f)  =  -H(f+R)  (1-4) 

so  that  the  unity  specification  about  f  =  0  automatically  accomodates  the  -1  specification 
about  f  =  R. 

Since  the  ideal  H(f)  behavior  is  unreachable,  one  may  define  a  frequency  limit  y  R/2 

and  a  magnitude  tolerance  6  which  require  N  and  b  choices  so  that 

n 

|  l-H(f)  |  ^  6  for  0  sf  s  y  R/2  (1-5) 

The  y  term  defines  a  passband  over  which  good  interpolation  is  achieved  and  the  6  term  con¬ 
trols  the  quality  level  of  the  interpolation.  The  passband  ripple  of  the  processor  (usually  not 
a  problem)  is  controlled  by  6  as  is  the  stopband  attenuation  (usually  controlling)  about  f  =  R. 

As  expected,  interpolator  complexity  (N  value  of  Equation  (1-3)  increases  with  increasing  y 
and  decreasing  6  .  in  later  portions  of  this  paper  interpolation  processors  will  be  used 
(  1,  Section  VI]  which  will  be  referenced  by  a  type  value  NT.  Type  selection  implies  a 
particular  (y  ,6  )  pair  where  the  6  value  is  usually  indicated  in  terms  of  stopband  attenuation 
level.  An  interestingclass  of  digital  filters  may  be  obtained]  1,  Section  IV  ]  from  the  im¬ 
pulse  response  of  the  interpolation  cascade.  (H(f)  of  Equation  (1-3)  has  an  elliptic  filter 
behavior  [  2,  Section  4.9J  .) 

The  well-known  duals  between  discrete  time  series  and  line  arrays  may  be  documented 
for  convenient  reference  with  the  aid  of  Figure  1-2.  The  frequency-domain  description  of 
the  discrete  time  series  of  Figure  l-2a  may  be  obtained  from  the  Discrete  Fourier  Transform 

-j27rk|- 

x<f)  =  2^  xk  e  (1_6) 

k 

where 

X(f+R)  =  X(f)  (1-7) 

so  that  quantization  (sampling)  in  the  time  domain  leads  to  a  periodicity  in  the  frequency 
domain. 
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If  the  symbol  T  is  also  used  to  represent  element  spacing  for  the  line  array  of  Figure 
l-2b,  the  far-field  pattern  Y(8  )  will  be  given  by 


Y(0) 


+  j  k 
e 


sine 

<W) 


(1-8) 


In  a  sense  (A/T)  represents  a  sampling  rate  since  this  ratio  equals  the  number  of  array 
elements  per  wavelength.  So  if  for  Equation  (1-8)  we  define 


f  =  sin  0 

(1-9) 

and 

R  =  A/T 

(1-10) 

We  obtain 

_  +j2?rk  R 

Y(t>  '  £  Vke 
k 

(1-11) 

where 

Y  (f  +  R)  Y(f) 

(1-12) 

which  reveals  that  sampling  in  the  aperture  domain  leads  to  periodicity  in  the  far  field 
domain  (sinO  space). 

While  there  is  only  a  trivial  sign  difference  between  Equations  (1-6)  and  1-11),  there 
are  physical  differences  in  the  situations  of  origin  which  should  be  kept  in  mind.  In  the 
array  case  f  values  outside  the  ±1  range  do  not  represent  physical  pointing  angles,  of  course, 
even  though  all  f  values  are  mathematically  acceptable.  The  selection  of  element  spacing 
(R  determination)  is  a  complicated  function  of  many  engineering  considerations  such  as 
mechanical  design,  maximum  scan  angle,  element  mutual  coupling  effects,  etc.  These 
factors  vary  considerably  from  one  application  to  another.  It  will  be  assumed  here  that  R 
for  the  line  array  is  pre specified  by  the  application  and  must  be  maintained  in  the  final 
design.  It  should  be  noted  from  Equation  (1-11)  that  the  far-field  pattern  in  sin  6  space  is 
directly  scaled  by  varying  the  element  spacing. 

In  the  next  section  interpolation  techniques  are  used  to  restore  the  sampling  rate  to  the 
prespecified  R  value  after  a  previous  operation  doubled  the  spacing  in  order  to  reduce  the 
sin  6  beamwidth  by  a  factor  of  two.  Thus  interpolation  is  used  to  fill  voids  in  a  stretched 
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array.  This  procedure  achieves  desired  beamwidth  and  sidelobe  level  specifications  by  a 
careful  insertion  of  new  elements  into  a  starting  "seed"  array  which  is  repeatedly  stretched 
and  filled.  Since  R  is  maintained,  the  aperture  size  is  expanded  by  the  process. 
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SECTION  II 


INTERPOLATION  CONCEPT 

The  pattern  synthesis  technique  to  be  described  uses  mathematical  procedures  which 
are  based  on  interpolation  methods  previously  discussed.  While  the  mathematics  and  the 
responsive  algorithms  may  be  described  in  a  unique  fashion,  there  are  two  physical  inter¬ 
pretations  of  the  process  which  justify  separate  expositions.  The  interpolation  concept 
follows  the  philosophy  of  the  previous  section  directly,  and  will  be  presented  first.  The 
"Auxiliary  Array"  concept  discussed  later  involves  a  much  different  interpretation  which 
represents  a  useful  alternate  physical  justification  of  the  design  procedure. 

Let  a  beamwidth  (unless  otherwise  stated  -3  dB  values  are  implied)  and  sidelobe  peak 
level  be  specified  for  an  array  for  a  given  element  spacing  or  R  value.  We  select  a  starting 
or  "seed"  array  which  has  the  requisite  sidelobe  level  but  (generally)  has  a  much  larger 
beamwidth  than  that  desired.  The  choice  of  seed  is  not  critical,  it  is  only  required  that  a  con¬ 
tiguous  sidelobe  region  exist  whose  peak  level  is  no  larger  than  that  specilied  for  the  final 
array.  A  convenient  seed  for  many  applications  may  be  represented  by  a  Chebyshev  design 
as  discussed  by  Dolph  and  others  (  3,4  ]  . 

The  seed  is  shown  in  Figure  2-l(a)  with  the  (prespecified)  element  spacing  T  which 
yields  the  required  R  value  (A/T).  The  element  weights  are  indicated  by  the  s^  terms.  The 
pattern  for  the  seed  is  sketched  in  Figure  2-l(c)  having  a  sidelobe  region  which  begins  at 
yR/2  in  sin  0  (f)  space.  The  seed  pattern  is  periodic  with  periods  equal  to  the  sampling 
rate  R  as  indicated  in  the  ST(f)  sketch.  We  now  stretch  the  seed  array  spacing  to  2T  as  in¬ 
dicated  by  the  s^  element  terms  of  Figure  2-l(b).  The  far-field  pattern  of  the  stretched 
seed  array  will  be  scaled  2:1  in  f-space  with  a  new  sampling  rate  of  (R/2)  as  indicated  by  the 
S2T(f)  plot  in  Figure  2-l(d).  Comparing  S2T(f)  with  ST(f)  in  Figure  2-1  it  is  seen  that  from 
an  absolute  f-scale  viewpoint,  the  stretching  has  narrowed  the  main  beam  which  is  desired, 
but  an  "echelon"  lobe  has  been  introduced  at  f  =  R/2  which  cannot  be  tolerated.  The  problem, 
therefore,  is  one  of  using  the  interpolation  processor  portion  of  Figure  1-1  to  produce  a  new 
set  of  weights  r^.  These  weights  will  be  placed  on  2T  centers  and  interleaved  with  the  s^ 
set  as  shown  in  Figure  2-l(b).  By  this  aperture  interpolation  process  we  hope  to  eliminate 
the  "echelon  lobe"  at  f  =  R/2  in  the  S2T(f)  function,  and  restore  the  aperture  sampling  rate 
to  the  presperified  R-value.  Note  that  the  interpolation  sequence  has  the  same  element 
spacing  as  the  "input"  or  "source"  array,  but  is  introduced  with  an  aperture  position  offset 
of  half  the  "input"  element  spacing,  (2T)/2  =  T,  The  comparisons  with  the  time  series  in¬ 
terpolation  (sampling  rate  multiplication)  process  of  Section  I  are  easy  to  make. 
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a)  Original  "Seed"  Array 


(v/2)  (R/2)  f\  (2-7)Rj 

1/2  (R/2)\___ 


(R/2) 


e)  H(2f),  [  1+H(2f)]  Behavior 


f)  Final  Output  Pattern  From  First  Processing  Level 


Figure  2-1.  Interpolation  Concept  Sketches  (Sheet  2  of  2) 


Consider  now  a  method  for  generating  the  weights  from  the  weights.  The  far- 
field  pattern  of  the  original  seed  is  given  by  Figure  2- 1(a). 
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where,  as  before 


f  =  sine  (2-2) 

R  =  A/T  .  (2-3) 

The  pattern  for  the  stretched  seed  is 
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which  simply  says  that 

S2T(f)  =  ST(2f)  .  (2-5) 

It  follows  at  once  from  Equation  (2-4)  and  Figure  2-l(b)  that  the  far -field  pattern  of  the 
array  is  given  by 


E2T(«.  £Tketi2'(k  +  1/2,<^)-  <2'6) 

k 

If  the  input  seed  values  s^  are  entered  sequentially  into  the  interpolation  processor  of 
Figure  1-1  they  will  appear  in  the  delay  element  in  reverse  order  with  the  lowest  index  values 
to  the  right  (convolution  effect).  Keeping  this  and  the  indexing  used  in  Figure  2-l(b)  in  mind, 
the  interpolation  weights  may  be  obtained  from 
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Since  a  symmetric  processor  is  to  be  used 
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allows  Equation  (2-7)  to  be  simplified  as 
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Substitute  Equation  (2-9)  into  Equation  (2-6) 
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Now  interchange  the  order  of  summation  to  get 
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Replace  k  by  k-n  in  the  first  interior  sum  and  by  k+n-1  in  the  second  interior  sum  to  yield: 
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Substitution  of  Equation  (2-4)  simplifies  to 
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E2T(f) 
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where 
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as  before.  Equation  (1-3)  permits  further  simplification  to 

E2T(f)  =  S2T(f)  H(2f)  -  (2-15) 

When  the  stretched  seed  and  the  interpolated  array  are  merged  (Figure  2-l(b)),  the  resultant 
pattern  is 

S^f)  S2T(f)  +  E2T(f)  =  S2T(f)[  1+H(2f)j  .  (2-16) 

From  the  discussion  in  Section  I  the  behavior  of  H(2f)  and  [  l+H(2f)J  is  sketched  in 
Figure  2-l(e).  H(2f)  remains  near  unity  out  to  y  R/4  then  follows  a  transition  region  from 
y  R/4  to  (2-y)  R/4.  At  this  f  value  the  -1  region  about  f  =  R/2  begins. 

The  (  1+H(2f)]  factor  starts  at  a  level  of  two.  At  the  center  of  the  transition  region 
(R/4)  the  value  is  one  (6-dB  down  from  peak  response).  The  "stopband"  region  begins  at 
(2-y)  (R/4)  with  a  peak  voltage  gain  of  6  .  Note  that  the  direct  6  suppression  value  in 
decibels  must  be  increased  by  6  dB  to  account  for  the  voltage  gain  of  two  in  the  "pass" 
region  of  (  l+H(2f)J  . 

The  final  output  Car-field  pattern  of  the  merged  array  is  obtained  by  taking  the 
product  S2T(f)  of  Figure  2-l(d)  with  the  J  l+H(2f)J  factor  of  Figure  2-l(e).  The  result  for 
Slp(f)  is  sketched  in  Figure  2-l(f).  The  echelon  lobe  of  S2T(f)  at  R/2  has  been  reduced  by 
the  process,  and  the  pattern  period  is  again  R  as  required.  The  new  sidelobe  regions  of 
SJj,(f)  may  be  identified  as  originating  (mainly)  from  the  seed  or  from  the  interpolation 
process.  These  regions  are  so  identified  in  Figure  2-l(f)  by  the  "s"  and  "I"  designations. 

In  Figure  2-l(f)  6  is  shown  at  a  lower  level  than  the  seed  sidelobe  level  for  delinea¬ 
tion  purposes  only.  Normally  these  two  sidelobe  levels  are  matched  in  design.  In  spite  of 
this,  sidelobe  peak  levels  only  are  directly  controlled,  with  sidelobe  peaks  normally  occur¬ 
ring  near  the  edge  of  the  main  beam  or  near  R/2.  Sidelobes  elsewhere  normally  drift  to 
lower  levels,  which  are  sometimes  brought  back  to  the  design  value  by  an  aperture  trunca¬ 
tion  process  to  be  described  later. 

One  can  now  logically  accept  the  merged  array  of  Figure  2-l(b)  as  a  new  "seed"  and 
repeat  the  whole  process  to  obtain  further  mainlobe  frequency  scaling.  The  basic  concept 
can  be  extended  to  any  number  of  levels.  In  fact,  the  interpolator  design  becomes  simpler 
at  the  later  levels.  Note  in  Figure  2-l(f),  for  example,  that  the  y  value  operative  here  is 
ono-half  that  encountered  in  the  true  seed.  Hence,  it  follows  from  the  discussion  of 
Figure  1-1  that  interpolator  design  becomes  simpler  at  each  succeeding  level  because  the 
design  y  is  cut  in  half  even  though  the  <5  specification  is  maintained. 
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(The  exceptionto  this  occurs  for  those  cases  in  which  the  main  beam  extends  beyond 
the  interpolator  yR/2  at  the  seed  level.  This  will  often  result  in  an  unacceptable  sidelobe 
peak  to  the  left  of  (2-y)  R/4  in  Figure  2-l(f).  A  second  level  y  larger  than  the  normal  half¬ 
value  will  be  needed  to  clear  this  up  (  1,  Section  3. 2j  .  After  this  special  "extended-gamma" 
second  level,  the  usual  half-gamma  sequence  may  be  followed  in  the  remainder  of  the 
cascade. ) 

In  the  multilevel  aperture  synthesis  procedure  the  output  at  each  level  contains  all  of 
the  elements  of  the  preceding  level  in  addition  to  the  new  interpolated  elements  generated 
at  that  level.  The  new  elements  are  generated  entirely  from  excitation  values  that  inhabit 
the  previous  level.  There  is  never  any  "reaching  back"  in  the  drive  matrix  beyond  the 
immediate  level.  An  element  once  created  at  a  level  appears  at  all  succeeding  levels,  in¬ 
cluding  the  final  one.  This,  if  one  looks  carefully  at  a  multilevel  synthesis  result,  the 
original  seed  will  always  be  found.  Note  also  that  the  operation  which  processes  2N  exci¬ 
tation  values  from  a  previous  level  (Equation  2-9)  to  produce  an  interpolated  value  for  the 
next  level  is  invariant  with  element  location  at  that  level.  The  design  changes  from 
level  to  level,  but  the  same  design  is  used  at  all  element  positions  within  a  given  level. 

Consider  now  an  alternate  use  of  the  interpolation  concept.  Let  the  seed  array  of 
Figure  2-2(a)  be  of  some  desired  ultimate  design  which  has  the  appropriate  beamwidth, 
sidelobe  levels,  and  element  spacing.  Now  let  it  be  required  that  a  second  array,  physically 
interleaved  with  the  first,  and  having  the  same  element  spacing,  T,  as  the  seed  array  be 
created  so  that  its  far-field  pattern  matches  that  of  the  seed  within  a  specified  tolerance 
over  a  specified  region  of  the  far  field.  (A  practical  need  for  this  "array  substitution" 
function  was  recently  brought  to  the  writer’s  attention. ) 

The  problem  now  becomes  one  of  deriving  a  set  of  element  weights  from  the  given 
set  s^  so  that  the  far-field  pattern  of  the  array,  E^.(f),  matches  the  far-field  pattern, 

S_(f)  of  the  s.  array.  If  the  logic  which  led  from  Equation  (2-1)  to  Equation  (2-15)  is 

1  K 

followed  without  the  doubling  of  T,  one  obtains  in  a  routine  manner  the  result 

ET(f)  =  H(f)  ST(f)  (2-17) 

where  H(f)  is  defined  by  Equation  (1-3)  as  before  and  Equations  (2-7)  and  (2-9)  are  valid 
without  alteration. 

If  the  r  currents  are  derived  from  the  s^  current  set  via  the  matrix  represented  by 

Equation  (2-9),  the  far-field  of  the  array  will  be  identical  to  that  of  the  s^  field  except 
for  a  possible  level  error  as  indicated  in  Figure  2-2(b).  The  gain  error  remains  small  out 
to  a  "cut  off  far-field  angle  which  correponds  to  y  R/2.  (This  is  essentially  a  limit  on  the 


2-7 


amount  of  beam  steering  that  may  be  done. )  The  closeness  of  the  gain  match  in  the  operat¬ 
ing  region  (6 )  is  a  function  of  interpolator  choice  as  is  the  y  value. 


h-  T  -H 
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(a)  Interleaved  Array  Geometry 
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(b)  Interleaved  Array  (t^  set)  Far-Field 
Response  Function 


Figure  2-2.  Displaced  Array  Substitution 
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Determination  of  the  weights  may  be  done  very  easily  using  the  techniques  of  this 
paper.  One  merely  inputs  the  s^  sequence  into  the  appropriately  selected  interpolation  sub¬ 
routine  as  described  in  Section  IV  and  the  interleaved  (s^,  t  )  set  of  Figure  2-2(a)  will  auto¬ 
matically  be  generated  by  the  program.  Additional  T-elements  with  intervening  zero-weight 
s-elements  will  automatically  be  added  at  the  array  ends  as  required,  to  align  phase  centers 
and  produce  the  desired  H(f)  performance.  The  number  of  such  edge  elements  is  kept  to  a 
minimum  by  this  design  procedure.  That  is,  for  a  given  y  and  <5  specification,  all  other 
design  procedures  will  require  a  larger  number  of  T-elements.  The  H(f)  function  character¬ 
izes  elliptic  filters  which  display  Chebyshev  behavior  in  both  the  pass  and  stopbands  [  2, 
Section  4.9]  .  Array  substitution  procedures  will  also  be  demonstrated  in  Section  VI  in  a 
different  context. 

Attention  is  turned  in  the  next  section  to  a  different  (but  perhaps  more  conventional) 
interpretation  of  the  "stretch-and-fill"  interpolation  operation. 
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SECTION  III 


AUXILIARY  ARRAY  CONCEPT 

In  the  interpolation  discussion  of  Section  II  and  Figure  2-1,  the  contributions  of  2N 
stretched-seed  weights  toward  the  generation  of  one  interpolation  weight  was  considered  in 
Equation  (2-7).  Here  the  process  is  reversed  and  it  is  noted  how  each  particular  seed  weight 
contributes  to  2N  different  interpolation  weights.  To  do  this  Equation  (2-7)  is  written  as  two 
separate  summations,  and  the  demand  is  made  that  the  s-subscript  remain  constant  at  some 
fixed  value,  say,  k.  For  the  first  sum  of  Equation  (2—7) 

k  =  k  -n  (3-1) 

and 


(3-2) 


(3-3) 

t  =  a  s  ,  n  =  1,  2,  —  N.  (3-4) 

k  +  n-1  k 

The  results  of  Equation  (3-2)  and  (3-4)  are  shown  in  Figure  3-1.  Each  seed  element 
of  Figure  2-l(b)  drives  an  auxiliary  array  of  2N  elements  with  the  excitation  matrix  shown. 
Nearby  seed  elements  similarly  generate  auxiliary  arrays  many  of  whose  elements  will 
coalesce  with  auxiliaries  of  other  seeds.  Because  of  superposition  we  may  still  treat  these 
auxiliaries  separately  in  order  to  demonstrate  an  underlying  principle  at  work. 

The  pattern  of  the  auxiliaries  will  first  be  found,  then  the  seed  element  will  be  brought 
down  into  place  (as  indicated  by  the  arrow  in  Figure  3-1)  and  the  overall  pattern  of  auxiliary 
with  seed  will  be  obtained.  From  Figure  3-1,  the  auxiliary  array  far-field  pattern,  A  (6), 
is  given  by  (normalized  to  unit  s_  current) 


r  =  a  s  ,  n  =  1,  2,  —  N . 

.-  -n  r 

k  -n  k 

For  the  second  sum,  the  k  value  must  be 
k  -  k  +  n-1 
so  that 
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Figure  3-1.  Auxiliary-Array  Generated  By  Each  Seed  Element 
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Substitution  of  Equations  (3-2),  (3-3)  and  (2-14)  yields 
N 

A(f)  =  £  cos  it  (2  n-1)  ^  . 

n=l 

Equation  (1-3)  may  be  used  to  show 
A(f)  =  H(2f) 


(3-5) 


(3-6) 


(3-7) 
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and  when  the  s_  element  is  brought  down  to  the  center  of  the  auxiliary  array,  the  effective 
k 

"element  pattern”,  E(f),  of  the  seed  and  its  directly  driven  auxiliaries  becomes 

E(f)  =  1  +  H  (2f)  .  (3-8) 

S  (f)  of  Figure  2-l(d)  may  be  thought  of  as  the  array  factor  for  the  stretched  seed. 
The  factor  [  1+H  (2f)  j  of  Equation  (3-8),  which  is  sketched  in  Figure  2-l(e),  becomes  the 
"element  pattern"  of  each  seed,  and  its  directly  associated  auxiliaries.  The  overall  far- 
field  pattern  is  the  product  of  the  array  factor  and  the  element  pattern  which  becomes  the 
S^,(f)  function  of  Figure  2-l(f). 

The  auxiliary  array  enhances  the  pattern  of  its  seed  in  a  nearly  uniform  manner  within 
the  main  beam.  In  the  echelon  lobe  region  about  f  =  R/2  (Figure  2-l(d)).  The  auxiliary  field 
cancels  the  seed  field  in  a  nearly  uniform  manner  in  order  to  produce  an  extended  null  in 
this  region. 

It  has  thus  been  demonstrated  that  the  interpolation  process  of  Figure  2-1  can  be 
decomposed  into  a  set  of  auxiliary  arrays,  one  for  each  seed.  The  reinforcement  in  the 
mainlobe  region  and  cancellation  in  the  echelon  lobe  region  has  been  explained  by  an  element 
pattern  anology.  This  viewpoint  forms  a  useful  alternate  interpretation  of  the  interpolation 
process.  A  further  development  of  the  auxiliary  array  concept  may  be  found  in  Appendix  B. 
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SECTION  IV 


PATTERN  SYNTHESIS  EXAMPLES 


Far-field  pattern  specifications  will  be  assumed  to  take  the  form  of  a  (-3  dB)  beam- 
width  and  a  peak  sidelobe  level.  The  synthesis  procedure  to  be  followed  involves  selection 
of  a  seed  array  having  the  requisite  sidelobe  suppression,  but  (normally)  a  much  larger 
beamwidth  than  the  specification  value.  An  interpolation  processor  of  the  type  described  in 
Figure  1-1  (or  the  equivalent  form  of  Figure  3-1)  is  selected  having  a  y  value  consistent 
with  the  main  beam  extent  of  the  seed  (as  per  Figure  2-1  and  associated  discussion).  The 
stopband  suppression  level  of  the  interpolator  is  chosen  to  meet  the  peak  sidelobe  specifica¬ 
tion  for  the  final  array,  and  is  usually  a  match  to  the  seed-sidelobe  level.  The  array  ele¬ 
ment  spacings  are  assumed  held  throughout  to  some  prespecified  value. 


The  desired  far-field  pattern  is  achieved  by  scaling  the  seed  pattern  repeatedly 
through  interpolation  stages  until  the  specified  beamwidth  is  achieved.  The  required  scaling 
factor  Fr  is  given  by 
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sin  B 
sin  B 


s/2 

r/2 


(4-1) 


where  B  is  the  seed  -3  dB  beamwidth  and  B  is  the  required  (specification)  -3  dB  beam- 
s  r 

width.  If  F  is  a  power  of  2,  then  log9  (F  )  will  yield  the  number  of  levels  of  interpolation 

that  will  be  required.  If  such  is  not  the  case  one  may  use  a  number  of  interpolation  levels 

in  excess  of  the  log  (F  )  value  and  then  one  may  resample  the  resulting  array  (select  every 
th  “  r 

J  element)  such  that 


F 


r 
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where  L  is  the  number  of  interpolation  levels  used  [  1,  Section  5.4]  .  An  example  of  such 
a  "J-derived"  array  will  be  shown  later. 

An  alternate  solution  to  this  problem  involves  variation  of  seed  design  so  that  the  re¬ 
sulting  value  of  Bg  in  Equation  (4-1)  will  yield  a  power  of  two  for  F  .  The  quantization 
forced  on  the  designer  (e.g.  number  of  seed  elements,  number  of  interpolator  coefficient 
pairs)  can  be  annoying,  but  such  is  the  nature  of  the  problem.  By  maintaining  a  discrete 
design  approach  throughout,  one  avoids  surprises  that  sometimes  occur  in  continuous/discrete 
conversions.  This  is  especially  helpful  when  the  final  array  contains  a  small  number  of 
elements. 
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The  synthesis  examples  to  be  given  here  do  not  necessarily  represent  optimum  designs. 
Readily  available  procedures  were  selected  for  convenience.  For  example,  the  interpolation 
processor  used  in  these  examples  was  originally  designed  for  other  purposes  [  1,  Section  4.2]  . 
As  shown  in  Table  4-1  four  interpolator  designs  were  used  having  nominal  suppression  levels 
of  -25,  -40,  -50,  and  -65  dB,  all  at  a  y  value  of  0.  8.  The  Chebyshev  seeds  were  then  de¬ 
signed  for  a  matching  mainlobe  region  (sidelobe  level  reached  at  f  =  yR/2)  and  sidelobe 
levels  as  close  to  the  interpolator  values  as  could  be  obtained.  As  seen  from  the  table,  the 
sidelobe  match  is  not  really  very  good  except  for  the  -25  dB  case.  It  might  also  be  men¬ 
tioned  that  the  interpolators  use  an  "extended  gamma"  design  at  second  level  [  1,  Section  3. 2] 
which  is  really  not  necessary  for  this  application.  Hence  the  interpolators  are  actually  some¬ 
what  overdesigned  for  the  intended  purpose. 

Some  examples  for  a  -25  dB  nominal  sidelobe  level  are  shown  first.  The  three- 
element  seed  pattern  which  is  used  for  this  set  is  shown  in  Figure  4-1.  In  all  of  these  plots 
a  normalized  sine  variable,  x=sin e/(R/2),  is  the  abscissa  variable.  The  sidelobe  region 
starts  at  x=0.8,and  has  a  null  at  x=0.86.  The  sidelobe  response  rises  to  the  specified  peak 
value  at  x=l.  The  x=l  value  is  a  "foldover"  value  for  the  pattern;  that  is,  the  patterns  are 
symmetric  about  x=l.  Thus  if  Figure  4-1  were  continued  to  x=2,  there  would  be  a  null  at 
x— 1. 14  and  the  mainlobe  region  centered  at  x=2  would  begin  x=l.  2.  If  the  seed  array  spacing 
were  doubled,  all  these  x-values  would,  of  course,  be  halved. 

The  pattern  resulting  from  a  single  interplation  level  is  shown  in  Figure  4-2.  This 
figure  will  be  examined  so  as  to  identify  first  the  effects  of  stretching  and  then  the  effects 
of  interpolation.  If  stretching  only  were  involved,  the  original  sidelobe  region  (0,8^x  ^  1) 
would  now  appear  from  (0.4  ^  x  ^  0.  5)  with  the  original  sidelobe  null  at  0.  86  appearing  at 
0.43.  These  effects  are  clearly  indicated  in  Figure  4-2.  The  sidelobe  continuation  in  x 
would  now  fall  at  (0. 5  -  x  ^0. 6)  with  a  null  at  x=G.  57.  The  pattern  would  then  build  from 
x  =0.6  to  a  second  mainlobe  peak  at  (0-dB  level)  at  x=l. 

The  effects  of  the  interpolator  may  now  be  understood  relative  to  Figure  4-2.  A  flat 
response  is  maintained  from  0  to  0.4.  A  transition  region  (see  Figure  2-l(e))  extends  from 
0.4  to  0.6  and  the  "stopband"  region  extends  from  0.6  to  1. 

The  nulls  in  the  pattern  beyond  x=0. 6  are  due  to  the  interpolator  as  is  the  pattern 
level  at  x=l.  (Without  the  interpolator  the  level  at  x=l  would  be  0  dB).  In  the  transition 
region  the  interpolator  contributes  0-dB  loss  at  x=0.4,  -6dB  loss  at  x=0.5,  and  its  design 
stopband  loss  (about  -26  dB  in  this  case)  in  the  region  (0.6  sx  £1).  This  behavior  can  be 
verified  by  close  examination  of  Figures  4-1  and  4-2. 
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TABLE  4-1 


PARAMETERS  USED  IN  EXAMPLES 


Figure 

No. 

No.  of 
Seed 

Elements 

Inter¬ 

polator 

Type 

Inter¬ 

polation 

Level 

No.  of 
Elements 

Element 
Count  After 
Truncation 

Beam- 

width 

(1) 

Side- 

lobe 

Spec 

Apper 
Effic  r| 

Comments 

4-2 

3 

4 

1 

15 

- 

17.38 

-25 

0.403 

4-3 

3 

4 

2 

35 

T- 

8.91 

-25 

0.339 

4-4 

3 

4 

3 

71 

- 

4.42 

-25 

0.336 

4-6 

4 

3 

1 

25 

- 

16.45 

-40 

0.270 

4-7 

4 

3 

2 

59 

- 

8.15 

-40 

0.230 

4-8 

4 

3 

3 

123 

- 

4.09 

-40 

0.221 

4-9 

1 

3 

1 

19 

- 

55.84 

-40 

0.114 

Omnidir. 

Seed 

4-10 

1 

3 

3 

99 

- 

13.46 

-40 

0.088 

Omnidir. 

Seed 

4-11 

5 

2 

3 

171 

- 

3.67 

-50 

0.177 

4-12 

6 

1 

3 

239 

~ 

3.27 

-65 

0.140 

5-3 

3 

4 

2 

35 

13 

10.03 

-25 

0.847 

Trunca¬ 
tion  of 

Fig  4-3 
Array 

5-4 

3 

4 

9 

4607 

13 

12.  25 

-25 

0.692 

J-De  rived 
Array, 

See  Text 

5-6 

4 

3 

3 

123 

47 

4. 14 

-40 

0.567 

Trunca¬ 
tion  of 

Fig  4-8 

5-8 

5 

2 

3 

171 

55 

3.64 

-50 

0.547 

Trunca¬ 
tion  of 

Fig  4-11 

5-9 

5 

2 

4 

347 

56 

3.67 

-50 

0.537 

Even  Ele¬ 
ment  Case, 
See  Text 

5-11 

6 

■ 

3 

239 

71 

3.27 

-65 

0.471 

Trunca¬ 
tion  of 

Fig  4-12 

a)  Far- Field  Pattern  Plot  Details 


»  ' 
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Notes: 


TABLE  4-1  (CONT) 
PARAMETERS  USED  IN  EXAMPLES 


NT 

y 

6  (dB) 

Design  Type 

1 

0.8 

-65 

Extended  y 

2 

0.  8 

-50 

Extended  y 

3 

0.8 

-40 

Extended  y 

4 

0.8 

-25 

Extended  y 

10 

0.6 

-50 

Normal 

b)  Interpolator  Data 


Figure 

No. 

No.  of 
Elements 

Beamwidth 

(1) 

Sidelobe 

Level 

(dB) 

3 

38.44 

-26.43 

4-5 

4 

33.52 

-42.07 

- 

5 

29.35 

-58.00 

6 

26.34 

-73.99 

c)  Seed  Pattern  Data 


1)  -3  dB  values  in  degrees  based  on  half  wavelength 
element  spacing. 


2)  Aperture  efficiency  =  rj 
No.  of  elements  =  K 

A 

Complex  element  weights 


K 

E 

k  =  1 


A 

w. 


K 

K  E 

k  =  1 


1 


The  pattern  after  a  second  interpolation  level  is  shown  in  Figure  4-3,  and  after  three 
interpolation  levels  is  demonstrated  in  Figure  4-4.  This  process  may  be  continued  indefin¬ 
itely  to  arrive  at  any  desired  final  beamwidth,  no  matter  how  small  this  may  be. 

It  will  be  noted  in  Figure  4-3  and  4-4  that  the  peak  sidelobe  levels  have  dropped  notice¬ 
ably  below  the  -25  dB  specification  value.  This  is  a  routine  occurrence  for  the  higher  inter¬ 
polation  cases  and  is  caused  by  two  main  factors.  First,  the  behavior  of  each  interpolator 
in  the  transition  region  adds  some  additional  attenuation  to  "specification  level"  sidelobes 
squeezed  into  that  region  by  the  stretching  process.  Secondly,  the  design  of  the  interpolators 
becomes  simpler  at  the  higher  interpolation  levels  since  the  effective  value  of  y  is  reduced  by 
a  factor  of  one-half  at  each  level.  This  generally  means  that  an  interpolator  having  a  smaller 
number  of  coefficient  pairs  (N  value  of  Equation  (2-7))  may  be  selected.  At  the  smaller  N  values 
quantization  becomes  a  factor.  An  interpolator  at  some  level  using  N  =  K  does  not  quite  pro¬ 
vide  enough  stopband  attenuation  while  the  N  =  K+l  design  provides  an  excess  of  echelon  lobe 
suppression.  One  always  opts,  of  course,  for  the  K+l  choice  in  designing  the  interpolator 
chain. 

The  use  of  an  overdesigned  interpolator  chain  produces  lower  than  specified  sidelobe 
levels  at  the  expense  (mainly)  of  the  creation  of  a  larger  than  necessary  number  of  "end- 
effect"  elements  at  the  array  edges.  Many  of  these  elements  have  such  small  excitation 
levels  that  they  may  be  removed  without  doing  serious  harm  to  the  pattern.  When  truncation 
is  properly  done,  the  sidelobe  peak  response  may  be  maintained  at  specification  level  while 
the  aperture  size  is  reduced  to  values  consistent  with  the  achieved  beamwidth  and  sidelobe 
level.  For  this  reason  the  element  counts  for  the  arrays  of  this  section  should  not  be  given 
much  significance.  This  topic  is  dicussed  in  more  detail  in  Section  V. 

Figure  4-5  shows  the  four-element  Chebyshev  pattern  used  as  the  seed  for  the  next 
set  of  examples  having  a  nominal  -40  dB  sidelobe  specification.  This  seed  yields  peak  side¬ 
lobe  levels  of  -42  dB  while  the  interpolator  has  a  -38  dB  stopband,  so  that  a  quantization- 
forced  mismatch  exists.  Figure  4-6  shows  the  pattern  after  one  level  of  interpolation.  The 
points  of  interest  here  may  be  examined  by  using  the  general  arguments  relative  to  Figures 
4-1  and  4-2  on  Figures  4-5  and  4-6.  The  more  complex  nature  of  the  interpolator  of  Figure 
4-6  is  indicated  by  the  number  of  nulls  between  x=0.6  and  x=l  as  compared  to  the  same 
region  of  Figure  4-2.  (This  is  the  echelon  lobe  region  in  both  cases. ) 

The  pattern  after  two  and  three  interpolation  levels  is  shown  in  Figures  4-7  and  4-8 
respectively.  These  last  two  figures  show  regions  of  very-deep  sidelobe  levels  in  addition 
to  a  peak  level  that  tends  to  fall  below  the  nominal  -40  dB  design  value. 
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Figure  4-4.  Pattern  of  3-Element  Seed  After  3-Interpolation  Levels  (-25  dB  Design) 


Figure  4-5.  Four-Element  Seed  Pattern  Used  For  -40  dB  Design 


Figure  4-6.  Pattern  of  4-Element  Seed  After  1 -Interpolation  Level  (  -40dB  Design) 


Pattern  of  4-Element  Seed  After  2-Interpolation  Levels  (-40  dB  Design) 


The  pattern  which  results  when  a  single -element  (omnidirectional)  seed  is  used  is 
shown  in  Figure  4-9.  The  seed  provides  a  "mainlobe"  region  at  all  angles,  hence  the  pattern 
at  first  level  reveals  the  workings  of  the  interpolator.  In  Figure  4-9  a  uniform  response  is 
maintained  from  x=0  to  x=0.4.  The  transition  region  extends  from  0.4  to  0.6,  while  the  stop- 
band  region  extends  from  x=0.6  to  x=l.  The  null  pattern  beyond  x=0.6  of  Figure  4-9  may  be 
compared  to  the  same  region  of  Figure  4-6.  The  stopband  behavior  of  the  interpolator  in  the 
echelon  lobe  region  of  Figure  4-6  is  identified  quite  easily. 

The  omnidirectional  seed  pattern  after  two  levels  of  interpolation  is  shown  in  Figure 
4-10.  The  "mainlobe"  region  is  scaled  down  by  a  factor  of  two  as  compared  to  Figure  4-9. 
While  these  last  two  figures  are  presented  as  oddities,  there  may  be  practical  cases  in  which 
a  uniform  response  in  the  main  beam  is  desired  accompanied  by  a  rapid  transition  to  the 
sidelobe  region.  Such  patterns  can  be  obtained  by  the  methods  outlined  here.  Additional 
comments  on  this  subject  will  be  made  in  Section  V. 

For  a  nominal  -50  dB  design  a  five-element  seed  is  used.  After  three  levels  of  inter¬ 
polation,  the  resulting  pattern  is  shown  in  Figure  4-11.  The  nominal  -65  dB  design  uses  a 
six -element  seed  whose  pattern  after  three  interpolation  levels  is  shown  in  Figure  4-12.  In 
both  cases  the  drop  in  peak  sidelobe  level  below  the  nominal  peak  value  is  to  be  noted.  The 
extended  areas  of  very  low  sidelobe  level  are  also  to  be  noted  in  these  two  figures.  These 
last  two  cases  are  presented  mainly  to  illustrate  some  points  to  be  detailed  in  Section  V. 

Table  4-1  summarizes  some  of  the  data  items  involved  in  the  pattern  figures  shown 
here  and  in  Section  V. 

Extensions  of  the  technique  to  greater  sidelobe  suppression  and/or  narrower  beam- 
width  are  routine.  The  few  cases  shown  here  are  meant  only  to  illustrate  design  principles 
which  will  be  further  defined  in  the  next  section. 
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Figure  4-9.  Pattern  of  Single-Element  (Omni)  Seed  After  1-Interpolation  Level  (-40  dB  Design) 


Pattern  of  Single  Element  (Omni)  Seed  After  2-Interpolation  Levels  (-40  dB  Design) 


SECTION  V 


CURRENT  DISTRIBUTIONS  AND  END-EFFECT  TRUNCATION 

In  previous  work  the  interpolation  processor  has  been  used  to  create  array  element 
weights  and  attention  has  been  focused  on  the  far-field  patterns  resulting  from  these  weights. 
In  this  section  attention  will  be  turned  to  the  nature  of  array  excitation. 

In  the  context  of  Figure  2-1  and  Equation  (2-7),  an  equivalent  interpolator  structure 
may  be  configured  as  shown  in  Figure  5-1.  This  arrangement  avoids  the  "reverse-order" 
entry  of  seed  values  as  required  by  the  Figure  1-1  arrangement.  It  is  convenient  to  think  of 
the  processor  of  Figure  5-1  as  being  moved  under  the  stretched-seed  array  producing,  at 
each  stop,  an  interpolation  value.  If  this  process  is  done  repeatedly  on  a  multilevel  basis, 
it  may  be  shown  that  the  number  of  coefficients  (elements)  at  level  L  is 

CL  =  (c0  ~  1>  2L  +  2  2(L+1"k)Nk-  2L+1  +3  (5-1) 

k=l 

where  CQ  is  the  number  of  elements  at  zero  level  (the  number  of  seed  elements),  and  N  is 
the  number  of  coefficient  pairs  at  level  k  in  the  interpolator  used. 


S.  .  S.  S.  8 .  _  S,  v. 

k-l  k  k+1  k+2  k*N 


Figure  5-1.  Interpolation  Processor  Realization 
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Since  every  weight  at  every  level  appears  in  the  final  array,  the  original  seeds  appear 
in  the  final  array.  We  define  as  "internal”  elements  of  the  final  array,  the  seeds  and  all 
elements  between  seed  elements.  By  "external"  elements  those  end  elements  are  meant 
which  are  found  outboard  of  the  outermost  seed  elements.  The  number,  Cj,  of  interior 
elements  is 

Cj  =  (C0-1)2L  +  1  (5-2) 

and  the  number,  C_  ,  of  exterior  elements  at  each  end  of  the  array  is 

IL 

CE  -  E  2<L+1'kV(2L-1)  (5-3, 

k=l 

so  that 

CL  =  Cj  +  2xC£  (5-4) 

In  typical  first-cut  designs  such  as  those  discussed  in  Section  V,  the  exterior  elements 
comprise  75  to  80  percent  of  the  array.  The  excitation  levels  in  these  elements  decrease 
rapidly  with  distance  from  array  center  and  a  high  percentage  of  these  elements  may  be 
discarded  in  an  aperture  truncation  operation  without  invalidating  the  far-field  pattern 
specifications. 

The  nature  of  this  aperture  illumination  ,  is  easily  explained  using  Figure  5-1.  The 
first  interpolator  coefficient  will  be  generated  at  level  one  when  the  a_N  weight  falls  under 
the  first  seed  element.  A  small  value  results  since  only  one  seed  contributes;  also,  the 
end  weights  of  the  interpolator  coefficient  set  tend  to  be  small.  (Note  too  that  alternate 
weights  in  this  region  are  zero  since  there  is  nc  seed  "brought-down"  between  interpolation 
values. )  Thus,  processor  "build-up"  and  "build-down"  regions  are  found  at  each  end  of 
the  level  one  array.  This  process  is  repeated  at  all  levels  except  that  beyond  level  one 
previous  "build-up"  regions  interact  first  with  the  processor  at  that  level,  and  the  tendency 
to  create  many  low-level  elements  at  the  array  ends  is  exacerbated. 

This  characteristic  of  the  procedure  is  also  heightened  by  use  of  overly  complex 
interpolation  processors  as  may  be  forced  on  the  designer  by  the  tyranny  of  quantization 
in  the  choice  of  N  as  discussed  in  the  previous  section.  An  important  advantage  of  this 
procedure,  however,  is  the  smooth  decay  of  the  illumination  function  that  always  results. 
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As  will  be  seen,  the  weight  function  is  not  quite  monotonic,  but  is  nearly  so  in  the  examples 
to  be  reviewed. 


By  nature  of  this  process,  the  weight  distribution  for  a  given  level  value  contains  the 
excitation  information  for  all  levels  less  than  the  one  shown.  (Selection  of  the  center  element 
and  every  other  element  at  level  L  yields  the  level  L-l  array,  etc. )  For  this  reason  aper¬ 
ture  functions  for  L=5  will  be  shown  even  though  far-field  patterns  for  lesser  L-values  are  of 
direct  concern. 

Figure  5-2  shows  the  level  5  weight  distribution  for  the  -25  dB  design.  In  these  figures 

bar  graphs  would  be  truly  representative  of  discrete -array  illumination.  However,  the  point- 

to-point  plots  presented  will  serve  since  the  point  density  is  high.  In  Figure  5-2  there  are 

th 

over  14  element  values  per  major  division.  Selecting  the  center  and  then  every  16u  element 
of  this  array  yields  the  pattern  of  Figure  4-2,  every  8th  element  gives  the  Figure  4-3  and 
selection  of  every  41*1  element  produces  the  far-field  pattern  of  Figure  4-4. 

The  level  two  pattern  of  Figure  4-3  is  generated  by  an  array  of  35  elements  taken  from 
the  weight  distribution  of  Figure  5-2.  By  symmetrically  truncating  this  35-element  array  so 
as  to  discard  low-weight  edge  elements,  a  13-element  array  having  the  pattern  shown  in 
Figure  5-3  results.  Truncation  is  very  easily  done  using  an  interactive  computer  program 
which  allows  pattern  effects  to  be  observed  as  elements  are  discarded.  One  merely  stops 
the  process  when  further  truncation  would  violate  a  specification. 

In  this  high-sidelobe  example  truncation  does  slightly  modify  main  beam  shape  but  the 
main  effect  is  in  the  sidelobes.  Comparing  Figures  5-2  and  4-3  shows  that  the  main  contri¬ 
bution  of  the  outboard,  low-level  elements  is  extra  sidelobe  suppression  as  argued  in  the 
previous  section.  Truncation  removed  22  of  35  elements  or  63%  of  the  reference  array. 

The  span  of  elements  kept  after  truncation  is  indicated  by  the  T-region  of  the  weight  plots. 

(In  all  of  the  weight  plots  shown  here  the  outermost  seeds  will  fall  at  about  0.4  and  0.6. ) 

The  opportunity  to  obtain  mainlobe  scaling  factors  which  are  other  than  a  power  of  two, 
was  outlined  previously,  and  is  demonstrated  in  Figure  5-4.  A  level  9  scaling  was  used 
(512:1)  which  was  later  resampled  by  a  J-step  of  163  to  obtain  an  effective  scaling  of 
512/163  =  3. 14  for  demonstration  purposes.  This  was  in  turn  truncated  to  13  elements 
yielding  the  far-field  pattern  shown  in  Figure  5-4. 

Figure  5-5  shows  the  aperture  illumination  for  the  -40  dB  design.  There  are  about 
25  elements  per  division  in  this  plot.  At  this  lower  sidelobe  level  a  relatively  smooth 
current  distribution  results.  Using  every  4th  element  from  this  aperture  gives  the  123- 
element  reference  array  of  Figure  4-8.  When  this  is  truncated  to  47  elements  the  pattern 
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of  Figure  5-6  results.  Comparing  with  Figure  4-8  again  reveals  that  truncation  restores 
sidelobe  levels  to  specification  values.  Keeping  47  out  of  123  elements  means  that  76  of 
123  or  62%  of  the  reference  array  elements  were  removed. 

Table  4-1  summarizes  some  of  the  data  items  involved  in  the  pattern  figures  of  this 
section. 

The  current  distribution  for  the  -50  dB  design  is  shown  in  Figure  5-7.  Selecting  every 
4**1  element  from  this  L=5  array  yields  the  171-element,  L=3  array  whose  far-field  pattern 
is  shown  in  Figure  4-11.  Truncation  of  this  reference  array  to  55  elements  produces  the 
pattern  shown  in  Figure  5-8.  Comparison  with  Figure  4-11  shows  a  dramatic  difference  in 
sidelobe  structure.  Deletion  of  116  elements  out  of  171  means  that  68%  of  the  reference 
array  elements  were  discarded. 

Routine  use  of  these  procedures  yields  an  array  having  an  odd  number  of  elements. 

If  an  even  number  of  elements  is  desired,  the  interpolation  concept  leads  quickly  to  a  simple 
procedure  for  odd-even  conversion.  One  simply  carries  the  design  to  an  interpolation  level 
one  step  beyond  that  originally  intended.  Then  a  2:1  censoring  is  done  such  that  the  central 
element  is  not  selected.  An  example  of  this  procedure  is  shown  in  Figure  5-9  which  may  be 
compared  directly  with  Figure  5-8. 

Figure  5-10  shows  the  current  distribution  for  the  -65  dB  design  which  contains  971 
elements  at  level  5,  The  pattern  of  the  level  three  array  derived  from  this  distribution  is 
shown  in  Figure  4-12.  When  this  array  is  truncated  from  239  elements  to  71,  the  pattern 
of  Figure  5-11  results.  Comparing  this  to  Figure  4-12  shows  that  even  though  the  peak  side¬ 
lobe  levels  have  been  raised  to  the  specification  value,  much  of  the  remaining  sidelobe 
region  is  suppressed  far  below  the  specification  value.  Retention  of  71  of  239  elements 
means  that  70%  of  the  reference  array  was  discarded  by  the  truncation  process. 

It  has  been  shown  that  this  synthesis  process  results  in  current  distributions  that  have 
desirable  tapers.  The  direct  (reference)  array  obtained  from  interpolation  contains  many 
elements  at  the  edges  with  small  weights.  Truncation  by  use  of  an  interactive  computer 
program  reduces  the  required  aperture  significantly  while  still  maintaining  prescribed  side¬ 
lobe  levels. 
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SECTION  VI 


STEERING  AND  FEED  MATRIX  ARCHITECTURE 

Up  to  this  point  the  problem  of  forming  a  beam  pattern  about  6  =0  having  prescribed 
properties  has  been  considered.  Attention  has  been  paid  to  the  number  of  elements  needed 
and  to  the  current  values  that  must  be  applied  to  each  element.  In  this  section  the  previous 
work  is  extended  to  encompass  the  ability  to  steer  the  beam  (shift  the  far-field  pattern  in 
sin  6  space)  and  means  for  creating  the  necessary  aperture  currents  which  produce  the 
steered  beam  pattern  are  developed.  A  unique  aperture  illumination  and  feed  matrix 
arrangement  is  described  which  accomplishes  the  beam  steering  function  without  the  use 
of  phase  shifters.  That  is,  the  feed  matrix  operations  are  entirely  real.  The  method  is 
presented  as  a  different  approach  to  beam  steering  which  gives  to  the  array  designer  an 
additional  option  from  which  to  choose. 

An  array  is  shown  in  Figure  6 -1(a)  having  an  odd  number  of  elements  with  weights 
"w  and  spacing  between  elements  equal  to  T.  For  the  moment  ignore  the  interleaved  array 
of  elements  which  also  has  element  spacing  T  and  whose  excitation  values  are  indicated  by 
gk<  Let  the"Wk  array  weights  be  derived,  for  example,  from  the  work  in  Section  V  which 
would  yield  a  symmetrically  driven  aperture  with  the  beam  center  at  6=0.  If  these  real 

A 

weights  are  replaced  by  the  complex  weights  W'k 

f0 

a  -j2*k  fg  f 

W,  =  W.  e  K  =  W.  cos  2^  k  —  -j  W,  sin  2ffk  —  (6-1) 
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The  far-field  pattern  becomes 
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(a)  Interleaved  Aperture 
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(b)  Quadrature  Current  Substitution  Details 


Figure  6-1.  Array  and  Feed  Matrix  (Sheet  1  of  2) 
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(c)  Feed  Matrix  and  Element  Drive  Modules 


(d)  Carrier  Signal  Generation 


Figure  6-1.  Array  and  Feed  Matrix  (Sheet  2  of  2) 


If  Eq  (f)  is  defined  as  the  far-field  pattern  existing  from  real  weights  (fQ  =0)  then 


E(f)  =  E0(f-f0)  (6-5) 

and  steering  has  been  accomplished. 

The  E(f)  pattern  is  real  and  is  generated  in  part  by  the  real  excitation  currents  and  in 
part  by  the  quadrature  or  imaginary  excitation  currents  of  Equation  (6-1).  The  corresponding 
components  of  the  far  field  are  given  by 


E  (f)  =  Wn  + 

R  0  k=l 


W. 


(f-'o) 

cos  2^rk  — — -  +  cos  27rk — — — • 
R  R 


and 


E,(f)  Lj  W. 

1  k  1 


(f-fo)  (f  +  fo)' 

COS  2  7T|<  -  -  COS  2^K  fj — 

R  K 


where 


E(f)  =  E  (f)  +  E  (f) 
R  1 


(6-6) 


(6-7) 


(6-8) 


The  real  currents  produce  a  pair  of  symmetrically  shifted,  real  beams  at  f  -  f  and 
f  +  fQ.  The  imaginary  currents  do  precisely  the  same  thing  except  for  a  polarity  reversal 
in  the  left-shifted  pattern  which  provides  for  the  cancellation  of  the  "image"  beam  and  re¬ 
inforces  the  beam  steered  to  f  =  fQ. 

Beam  steering  clearly  requires  the  effects  of  complex  currents  at  the  elements  and 
it  is  normal  practice  to  employ  phase  shifters  in  excitation  matrices  of  beamsteering 
phased  arrays.  We  propose  here  to  eliminate  the  need  for  phase  shifters  by  utilizing  a  pair 
of  interleaved  arrays,  each  fed  with  currents  of  one  phase  only.  The  set  of  Figure 
6 -1(a)  will  be  fed  entirely  with  real  or  0°  phase  currents,  while  the  gk  set  will  be  fed  with 
imaginary  or  90°  phase  currents.  Polarity  reversals  in  both  sets  due  to  negative  real 
weights  will  be  allowed.  The  currents  in  the  Wk  set  will  be,  by  our  option,  the  real  com¬ 
ponent  of  of  Equation  (6-1) 

f 

Wk  =  Wk  COS  2n  k  T  (6~9) 

The  imaginary  component  of  current  formerly  at  W. 

f  k 

qR  -J  Wk  sin  2*k  (6-10) 
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has  been  omitted.  An  attempt  will  be  made  to  excite  nearby  elements  so  that  the  far- 
field  will  be  unaware  of  the  transfer  of  quadrature  current  from  to  its  quadrature -excited 
neighbors.  This  cannot  be  done  exactly,  but  as  will  be  seen,  it  can  be  done  with  sufficient 
accuracy  in  a  practical  manner  if  proper  interpolation  procedures  are  followed. 

All  quadrature  currents  must  have  an  odd-symmetric  counterpart  in  the  array  since 
the  far-field  pattern  is  to  be  maintained  real.  A  unit  quadrature  current  (+j)  at  distance  x 
from  array  center  (along  with  its  counterpart)  produces  a  far-field  phasor 


f 

E.  =  -2  sin  2?rx  r 
J 

The  field  of  the  original  quadrature  current  of  Equation  (6-10)  is  then 

f 

E^  (f)  =  2  sin  2?rk  — sin  2ff  k  ^ 
which  becomes  the  target  of  the  approximation  process  which  follows. 


(6-11) 


(6-12) 


The  far-field  pattern  of  Equation  (6-12)  will  be  synthesized  by  N  pairs  of  symmetri¬ 
cally  excited  elements  (gk  -  gk_j),  (gk+1  =  gk_2),  — ,  (8k  ^.4  =  Sk_N>etc.  as  detailed  in 
Figure  6-l(b).  The  far-field  pattern  of  these  N  pairs  (and  their  odd-symmetric  counter¬ 
parts)  is  given  bv 


N 


\<f)  =  in4?  -4  gk+n-i cos  27r  1  j sin  27rk  I 


£  -4  cos  27r^=~^  --  f  sin  2  7rkt  (6-13) 

Now  it  is  desired  that  A.  (f)  of  Equation  (6-13)  agree  as  closely  as  possible  with  E  (f)  of 
Equation  (6-12),  an  identity  is  not  possible  so  we  write 


N  f 

-4Skfn_1cos27r^-  1 t 2  Wk  sin  2*k  -f 


(6-14) 


For  any  given  k  value  and  steering  angle  fQ,  the  right-hand  side  is  a  constant  (Wk  is 
also  prespecified).  For  illustrative  purposes  we  normalize  by  defining 
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then  Equation  (6-14)  may  be  written  in  the  familiar  form  of 


H(f)  =  E  bn  cos  2tt  |  »  1 


(6-16) 


which  is  H(f)  of  Equation  (1-3). 

A  set  of  coefficients  bR  can  be  found  which  will  hold  the  approximation  of  Equation 
(6-16)  to  within  a  specified  error  out  to  a  design  cutoff  of  f  =  yR/2  as  extensively  discussed 
in  previous  sections. 

Suppose  now  that  an  an,  bR  set  is  selected.  The  g-values  required  in  the  neighbor¬ 
hood  would  be  given  by,  (using  Equation  6-15) 


’k+n-l 


*n  (  ~Wk  sin  27rk  -g-  )  •  n  =  1, 2, - N 


(6-17) 


which  would  yield  (gk,  gkl),  (gk+1>  gk_2>»  —  (gk+N1>  gk  NJ)  because  of  the  assumed 
symmetry.  Thus  the  single  computed  value  (-Wk  sin  2?rk  f^/R )  contributes  to  2N  g-values 
through  the  an  weights.  Worked  the  other  way  around,  this  means  that  each  g-element  of 


the  array  will  be  fed  from  2N  of  the  qk  "missing  currents"  of  Equation  (6-10). 
expression  for  the  g-currents  becomes 


N  r 

E  \ 

n  =  1 


f0  f0 
~Wk+n  Sin  27r  (k+n)  ~R  "  Wk+l-n  Sin  27r(k+1“n>  R- 


Thus,  the 
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which,  not  surprisingly,  amounts  to  an  interpolation  process  over  a  neighboring  set  of 
"missing"  quadrature  currents  (see  Equation  (6-10)).  The  real  currents  are  available 
without  ceremony  from  the  simple  calculation  of  Equation  (6-9)  and  are  used  to  feed  the 
original  Wk  elements  as  described  before. 

It  is  interesting  to  note  that  this  process  of  converting  a  set  of  array  samples  requir¬ 
ing  complex  operations  for  their  generation  to  a  set  of  array  samples  permitting  real  oper¬ 
ations  in  their  generation  is  identical  to  a  signal  processing  procedure  treated  earlier  by  the 
writer  (  Ref.  1,  Sections  5. 1,  5.  2,  see  especially  Equation  (4-4))  .  In  fact  the  computer 
program  needed  to  implement  conveniently  the  feed  matrix  of  Figure  6-l(c)  is  available 
[  Ref  1,  Section  6.1,  subroutine  CORE,  complex-to-real  transformation)  . 

(In  signal  processing  work  sampled  data  is  usually  found  in  one  of  two  formats.  A  band¬ 
pass  signal  may  be  downshifted  to  a  low  (baseband)  frequency,  and  this  translated  version  is 
sampled  to  yield  "real  samples".  Alternately  a  cosine  and  sine  detection  near  band-center 
of  the  signal  may  be  done  and  the  resulting  lowpass  I  and  Q-components  may  then  be 
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sampled  to  yield  "complex  samples",  it  often  occurs  in  signal  processing  work  that  trans¬ 
formations  from  one  sampling  format  to  the  other  are  required.  It  was  shown  in  the  reference 
previously  cited  that  these  transformations  are  in  reality  interpolation  problems.  Subrou¬ 
tines  for  complex-to-real  (CORE)  and  real-to-complex  (RECO)  conversions  are  presented 
and  described  in  the  previous  reference. ) 

A  few  comments  now  on  the  synthesis  of  the  required  odd-symmetric  far  field  by  the 
g^  set  as  described  in  Equations  (6-13  through  6-16).  Equations  (6-6  and  6-7)  show  that 
failure  to  produce  E^(f)  will  leave  the  f  +  fQ  portion  of  E^(f)  uncancelled  and  an  "image"  beam 
will  appear.  Substitution  of  Equation  (6-17)  into  Equation  (6-13)  yields,  using  also  Equations 
(6-16  and  6-7), 

*0  f 

Ak(f)  11(f)  2  Wk  sin  2nk  sin  2nk  ^  (6-19) 

and  therefore 

£  v  £  wk 

k  k 

When  this  is  compared  to  Equations  (6-12)  and  (6-7)  it  will  be  seen  that  the  surrogate 

quadrature  currents  are  doing  exactly  the  right  thing  except  for  level  (amplitude)  errors. 

From  f  0  to  f  yR/2  H(f)  has  a  maximum  error  (deviation  from  unity)  of  ±  6  as  discussed 

in  Section  I.  For  f  >  yR/2,  H(f)  ceases  its  valiant  struggle  for  unity  and  quietly  goes  to  zero 

at  f  R/2.  The  region  beyond  f  =  yR/2  will  be  an  errored  region  in  which  the  ED(f)  of 

n 

Equation  (6-6)  will  prevail  undiminished  while  the  g^  set  will  only  be  able  to  supply  a  frac¬ 
tion  of  the  Ej(f)  of  Equation  (6-7).  Real  harm  will  occur  in  such  cases  only  if  the  beam  is 
steered  to  an  f  value  in  (yR/2,  R/2)  or  (-  yR/2,  -R/2).  The  true  main  beam  may  appear 
somewhat  distorted  in  the  process,  but  an  uncancelled  image  will  appear  at  the  other  f-space 
extremity.  The  failure  of  the  image  to  cancel  cannot  be  tolerated.  Steering  to  ±R/2  repre¬ 
sents  an  ultimate,  obviously.  Interplators  that  will  give  us  80%  to  90%  of  this  space  are 
readily  available.  The  main  price  paid  for  higher  y  is  (usually)  a  modest  extension  of 
aperture. 

In  the  working  region  (-yR/2  *  f  <  yR/2)  the  failure  of  the  g^  set  to  cancel  the  image 
main  beam  will,  in  essence,  tend  to  violate  the  sidelobe  specification  for  the  array.  The 
Figure  6-1  interpolator  must  therefore  have  a  stopband  specification  at  least  equal  to  the 
sidelobe  level  desired. 

These  matters  are  demonstrated  in  Figures  6-2  through  6-8.  The  feed  matrix  of 
Figure  6 -1(c)  was  programmed  using  the  signal-processing  subroutine  CORE  (the  polarity 
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Reference  Pattern  Steered  to  +0. 2R.  Beam  Steering  Interpolator 


Reference  Pattern  Steered  to  +0.  2R.  Beam  Steering  Interpolator 


Reference  Pattern  Steered  to  +0. 35R.  Beam  Steering  Interpolator 


reversals  provided  by  CORE  [  Ref  1,  Appendix  A]  must  be  eliminated  for  this  use)  for  imple¬ 
mentation  of  the  interpolator  section  of  Figure  6-l(c)  as  previously  discussed.  The  reference 
pattern  chosen  for  this  group  of  examples  was  the  -50  dB  pattern  of  Figure  5-8.  This  55- 
element  pattern  was  chosen  because  the  low  sidelobe  levels  would  make  immediately  obvious 
any  performance  deviations  of  the  feed  matrix.  As  will  be  seen  the  structure  of  Figure  6-l(c) 
performs  exactly  as  advertised. 

Figure  6-2  shows  the  reference  pattern  of  Figure  5-8  in  the  plotting  format  of  this 
group.  The  beam  is  centered  in  the  broadside  direction  with  f  -  ±R/2  at  the  plot  extremities. 
Thus,  these  plots  show  an  entire  far-field  period.  The  -55  dB  sidelobe  at  f  =  ±R/2  should  be 
noted.  With  fQ  set  to  zero,  it  follows  from  Figure  6-l(c)  that  all  the  g^  currents  are  zero, 
while  the  W'^  currents  equal  the  original  W^  weights.  Thus  the  array  is  restored  exactly  to 
the  unsteered  format  analyzed  in  Section  V. 

Figure  6-3  shows  steering  for  fQ  =  -0.  2R.  The  -55  dB  sidelobe  has  moved  down  from 
f  0.5R  to  f  0.3R.  Note  the  "image"  of  the  main  lobe  at  f  -  +0.2R.  The  beam  steering 
interpolator  is  the  NT  =  2  choice  which  provides  a  stopband  attenuation  of  about  -50  dB.  This 
is  well  demonstrated  by  the  "sidelobe"  response  at  f  =  +0.  2R.  Comparing  these  last  two 
figures  shows  that  in  all  major  features  a  pattern  rotation  has  indeed  been  implemented  by 
the  feed  matrix. 

The  beamsteering  interpolator  y  of  0.  8  at  NT  =  2  indicates  that  the  fQ  operating  region 
of  the  matrix  extends  to  ±0.8  x  0.5R  =  ±0.4R.  Figure  6-4  shows  steering  to  fQ  =  +0.35R. 

The  -55  dB  sidelobe  has  moved  from  -0.  5R  to  -0. 15R  and  the  "image"  of  the  main  beam  at 
f  -0. 35R  is  clearly  evident.  (The  shape  of  the  image  will  change  depending  upon  its  exact 
location  within  the  H(f)  pattern. )  The  pattern  shown  here  near  the  design  edge  of  the  matrix 
is  proper  and  correct. 

Figure  6-5  shows  steering  beyond  the  design  limits  to  fQ  =  +0.45R.  In  this  region  the 
surrogate  quadrature  currents  can  no  longer  produce  the  properly-sized,  far-field  contribu¬ 
tion.  The  main  lobe  loses  some  level  (not  noticeable  with  plot  normalization).  The  main 
effect  however  appears  in  the  image  region  in  which  two  relatively  large  pattern  components 
can  no  longer  produce  the  small  difference  which  was  held  to  a  level  below  -50  dB  in  the 
operating  region.  The  "sidelobe"  level  observed  here  is  of  the  order  of  -14  dB.  This  is  a 
rather  clear  indication  of  a  serious  deviation  from  the  norm. 

In  Figure  6-6  the  beam  steering  interpolator  is  changed  to  an  NT  =  4  design  which  has 
a  y  of  0.  8  but  only  a  -25  dB  stopband  specification.  The  "image"  beam  at  f  =  0. 2R  in  re¬ 
sponse  to  the  main  beam  at  fQ  =  +0.  2R  clearly  shows  the  design  mismatch.  With  the  excep¬ 
tion  of  this  rather  serious  lapse,  the  remainder  of  the  pattern  of  Figure  6-6  is  quite 
satisfactory. 
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Figures  6-7  and  6-8  use  a  beam  steering  interpolator  choice  of  NT  =  10  which  has  a 
-50  dB  stopband  design  but  a  yof  only  0.6.  Thus  the  operating  limits  here  on  fQ  are 
±0.6  x  .  5R  =  ±0. 3R.  Figure  6-7  shows  the  pattern  for  a  pointing  angle  of  fQ  =  +0. 2R.  Normal 
operation  results  with  the  image  beam  well  below  the  -50  dB  specification  level. 

Figure  6-8  shows  the  same  beamformer  steered  to  +0.35R,  which  is  beyond  the  ±0. 3R 
design  limit.  As  expected,  the  problem  area  is  at  f  =  -0.35R  where  a  -28  dB  image  beam 
appears. 

These  examples  confirm  the  design  philosophy  which  led  to  the  design  of  the  feed 
matrix  of  Figure  6-l(c).  The  extension  of  aperture  needed  to  accomodate  the  interleaved 
elements  has  been  of  the  order  of  20%  for  the  data  shown.  These  are  raw,  untruncated 
results  and  should  not  be  given  any  special  significance.  Quick  looks  at  the  element-weight 
listings  indicate  that  considerable  truncation  may  be  done  as  was  the  case  in  Section  V.  An 
estimate  of  the  aperture  extension  required  for  the  above  cases  might  be  5  to  10%.  This 
matter  will  be  more  thoroughly  studied  when  time  and  circumstances  permit. 

Extensions  of  this  work  to  the  case  of  the  planar  array  are  treated  briefly  in 
Appendix  A. 
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SECTION  VII 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

This  work  represents  a  different  approach  to  the  design  of  discrete  arrays.  The  sub¬ 
ject  area  has  been  intensively  pursued  over  the  years  with  outstanding  results  by  many 
capable  investigators.  The  main  value  of  this  particular  work  lies,  perhaps,  in  the  exposi¬ 
tion  of  new  points  of  view  on  an  old  problem.  Such  transients  in  established  thought  patterns 
can  quite  often  produce  beneficial  effects. 

The  data  shown  in  Table  4-1  indicates  that  the  design  procedure  is  generally  competi¬ 
tive  with  prior-art  approaches  in  terms  of  aperture  size  and  efficiency.  The  procedure  is 
based  on  discrete  arrays  from  the  outset.  This  eliminates  problems  that  can  sometimes 
arise  when  continuous  distributions  are  sampled  to  obtain  element  weights. 

The  procedure  also  lends  itself  to  implementation  using  interactive  computer  tech¬ 
niques.  This  is  especially  true  in  the  truncation  operation  described  in  Section  V.  In  many 
cases  the  reference  arrays  obtained  directly  from  the  interpolator  cascade  have  extended 
sidelobe  areas  of  very-low  level.  Truncation  tends  to  fill  these  regions  back  up  in  exchange 
for  a  smaller  aperture  and  higher  efficiency.  In  some  cases  a  larger  than  minimum  aper¬ 
ture  and  a  lower  than  maximum  efficiency  may  be  readily  accepted  in  exchange  for  extra 
sideloix'  suppression  in  extended  regions.  Interactive  computer  operation  usage  of  these 
procedures  can  be  most  helpful  in  tailoring  an  array  design  to  a  specific  practical  application. 

It  also  ap|x?ars  that  the  freedom  given  to  the  array  in  terms  of  sidelobe  variation  (only 
the  peak  value  is  directly  manipulated  by  the  design)  may  contribute  to  the  well-behaved 
nature  of  the  current  distributions  which  result.  It  is  the  nature  of  this  procedure  that 
current  distributions  will  always  taper  smoothly  to  zero  at  the  aperture  edges  regardless  of 
the  beamwidth  or  sidelobe  specification  imposed.  A  minor  advantage  of  the  technique  may 
lie  in  the  fact  that  weight  calculations  for  any  pattern  specification  are  straightforward 
almost  to  the  point  of  being  trivial. 

The  feed  matrix  and  aperture  combination  developed  in  Section  VI  has  the  immediate 
advantage  of  novelty.  Conventional  beam  steering  employs  no  amplitude  control,  but  provides 
individual  phase  adjustment  for  each  element.  The  system  of  Section  VI  uses  one  phase 
shifter  but  requires  level  (and  polarity)  control  at  each  element.  For  constant-envelope  wave¬ 
forms,  the  phase  shifter  approach  is  a  compatible  one.  If  more  general  types  of  signal  wave¬ 
forms  are  to  be  handled,  then  the  options  will  have  to  be  reexamined.  At  the  least,  the  beam¬ 
forming  system  of  Section  VI  gives  to  the  designer  one  more  option  from  which  to  make 
a  selection. 
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The  basic  concept  demonstrated  in  Section  VI  involves  the  use  of  interpolation  tech¬ 
niques  to  generate  currents  for  use  in  a  displaced  element  set  to  match  the  far-field  pattern 
of  an  original  element  set.  The  specific  use  made  of  this  concept  permitted  separation  of 
in-phase  and  quadrature  current  radiators  in  an  array.  However,  this  certainly  does  not 
exhaust  the  possible  uses  of  the  basic  concept  demonstrated. 
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APPENDIX  A 


EXTENSION  OF  ARRAY  SUBSTITUTION 
PROCEDURES  TO  TWO  DIMENSIONS 

The  array  substitution  procedures  of  Section  VI  will  be  extended  to  cover  the  planar  array 
case  in  the  following.  Figure  A-l(b)  shows  the  coordinate  system  chosen  here  for  the  discussion. 
The  far- field  point  E  lies  at  a  distance  P  from  the  origin  with  azimuth  angle  0  and  elevation 
angle  0 .  The  array  is  assumed  to  be  in  the  y-z  plane  with  element  spacing  T  as  shown  in  Figure 
A- 1(a).  Let  the  array  be  in  the  first  quadrant  of  the  y-z  plane  and  let  an  element  location  be 
specified  bv  its  (y,z)  coordinates. 

If  d  is  defined  as  the  distance  from  the  element  to  the  far-field  point,  E,  some  simple 
geometry  yields,  for  P  -*»  « 

p  -  d  =  y  cos  &  sin  0  +  z  sin  0  (A-l) 

The  phase  advance  due  to  distance  of  an  element  at  (y,  z)  relative  to  an  element  at  the 
origin  is 

O  =  — ( —  (y  cos  0  sin  0  +  z  sin0  )  (A-2) 

a  \ 


Consider  now  an  element  of  the  original  array  W.,  which,  for  convenience,  will  be  nor- 

malized  to  unit  current.  It  is  also  convenient  to  consider  the  coordinate  system  origin  to  be 

relocated  to  the  W..  location.  It  is  now  desired  that  the  far-field  pattern  of  the  unit  current 
Jk 

W  k  be  duplicated  as  closely  as  possible  by  a  substitution  element  set  g-k  in  the  W.^  neighbor¬ 
hood.  If  the  phase  of  the  substitution  element  set  far-field  pattern  is  required  to  match  the 
phase  of  the  W.^  radiation  at  all  far-field  points,  the  g.k  must  be  chosen  in  groups  of  four  with 
each  member  of  the  group  (quad)  driven  by  the  same  current;  that  is,  for  every  first  qua¬ 
drant  element  at  (y,z)  (relative  to  W.k),  one  needs  three  additional  elements  at  (-y,z),  (-y,-z), 
and  (v,-z)  all  driven  with  the  same  current,  I.k> 

With  y  and  z  measured  from  W^k,  the  far-field  pattern  of  any  such  quad  will  be 


E..  (  «,  0  )  =  21 

Jk 


jeos  [-^  <y 


cos  0  sin  0  +  z  sin0  ) 


+  cos  — r  (-y  cos  0  sin  9  +  z  sin 


•'ll 


(A-3) 


where  y  and  z  are  the  coordinates  of  the  first-quadrant  element  and  I^k  is  the  common  current 
drive  to  the  four  substitution  elements. 


O  *  Original  element  set 
■+•  *  Substitution  element  set 
R  »  T/X 


(b)  Coordinate  System 
Figure  A- 1.  Planar  Array  Description 


Define 


f  =  sin  0 

(A -4) 

g  A  sin  <f> 

(A-5) 

\=.2a 
r  r 

(A-6) 

and  let  (Figure  A- 1(a)) 

y  =  (2j-l)  | 

(A- 7) 

z=  (2k- 1)  | 

(A-8) 

and  require  that 


I.,  =  a.  a,  . 

Jk  J  k 

Equation  (A- 3)  may  now  be  rewritten  as 
E.r  (f,  g)  =  b.  bk  cos  |  -g-  (2j— 1 )  /l-g2  *' f  ] 

x  cos  [  ^  (2k- 1)  gj  . 


(A-9) 


(A-10) 


This  equation  gives  the  far-field  contribution  of  four  particular  elements  which  are  symmet¬ 
rically  located  about  W.,  ,  and  identified  by  the  row,  column  coordinates  of  the  first  quadrant 

1K 

member  of  the  quad.  If  N  rows  of  g^k  elements  above  and  below  W^k  are  to  be  used  in  the  sub¬ 
stitution  process  along  with  N  columns  to  the  right  and  left  of  W^k>  then  the  total  far-field  pat¬ 
tern  of  the  substitute  array  for  W.k  will  be 


N  N 


E(f,g)  =  la  la 


Vf-Sl 


j=  1  k  =  1 
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"  ^  b.  cos  I (2J-1)  l/7f|  cos  l-g-gk-Dgl 

J-l  ’ 


so  that 


E  (f , g)  =  H  (f  j/T?)H 


(g) 


(A-ll) 


where  H  (f)  is  the  ubiquitous  function  defined  in  Equation  (1-3) .  It  is  useful  to  consider  Equa¬ 
tion  (A-ll)  to  be  an  effective  element  pattern  of  the  substitution  element  set  whose  array  factor 
is  obviously  the  same  as  that  of  the  original  element  set.  It  then  follows  that  the  far-field  pat¬ 
tern  of  the  substitution  element  set  is  that  of  the  original  array  times  the  E  (f,  g)  function  of 
Equation  (A-ll). 

It  should  be  recalled  that  in  the  above  development  a  particular  excitation  procedure  was 
involved  for  a  set  of  substitution  elements  in  the  neighborhood  of  an  original  element  in  an  at¬ 
tempt  to  duplicate  the  far-field  function  of  that  original  element.  Equation  (A-ll)  shows  that 
the  substitute  far-field  pattern  matches  the  original  without  error  in  phase.  The  error  in 
level  is  controlled  by  the  H  function  product  of  Equation  (A-ll. )  In  the  x-y  plane,  0  =  0,  g=  0 


E  (f,  0)  =  H  (0)  H  (f) 


(A-12) 


In  the  x-z  plane,  0  =  0,  f  =  0  and 
E(0,g)  =  H  (o)  II  (g) 


(A- 13) 


Since  H  (0)  =  1,  the  behavior  in  these  principal  planes  is  similar  to  that  of  the  line  array 
as  discussed  in  Section  VI.  For  a  given  N  value  and  corresponding  an  set  of  Equation  (1-3), 
the  behavior  of  H(f)  is  known.  Therefore,  Equation  (A-ll)  may  be  used  to  determine  the  per¬ 
formance  level  of  the  substitution  array  for  any  f,  g  combination. 

Steering  of  the  planar  array  will  require  complex  element  currents.  If  desired,  the  real 
components  may  be  kept  in  the  original  array  and  the  complex  components  may  be  placed  in  the 
substitution  element  set  as  was  done  for  the  line  array  of  Section  VII  to  avoid  the  need  for  phase 
shifters.  Once  the  relationship  of  Equation  (A-ll)  is  established,  the  equivalences  between  the 
line  and  planar  array  are  so  close  as  to  make  elaboration  superfluous. 

Equation  (A-9)  was  based  on  a  unit  current  in  the  element  W.,  for  which  the  substitution 
was  being  done.  The  indices  j,  k  were  a  measure  of  the  y-and  z-distances  from  W.,  as  is 

JK 

clear  from  Equations  (A-7)  and  (A-8).  If  one  considers  the  total  current  in  g.,  resulting  from 
substitution  operations  targeted  at  a  group  of  original  elements  W4.  in  the  g.,  neighborhood, 
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there  results 
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r  j-r  +  1,  k-q  +  1  r  j  +  r,  k-q  +  1 


(A- 14) 


This  is  a  double  interpolation  process.  The  first  interior  sum  on  r  represents  an  interpo¬ 
lation  to  obtain  intermediate  element  currents  in  the  rows  above  g.,  (compare  Equations  (2-9) 

jk 

and  (A-14)).  The  second  interior  sum  on  r  obtains  the  intermediate  element  currents  in  the 
rows  below  g  .  The  outer  sum  over  q  in  Equation  (A-14)  represents  an  interpolation  operation 
on  the  intermediate  row  values  in  order  to  obtain  the  intermediate  column  values  which  are  the 

It  follows  that  Equation  (A-14)  defines  an  efficient  feed 

j*' 

matrix  algorithm. 


actual  substitution  element  currents  g^. 
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APPENDIX  B 


A  SWITCHABLE  END- FIRE  ELEMENT  FOR  USE  IN  THREE-DIMENSIONAL  ARRAYS 


A  particular  extension  of  the  work  of  Section  HI  is  presented  here.  The  array  of 
Figure  3-1  is  redrawn  in  Figure  B-l(a)  while  the  far-field  pattern  of  Figure  2-l(e)  is  shown 
for  convenience  in  Figure  B-l(b).  The  array  of  Figure  B-l(a)  has  a  pattern  which  is  very 
nearly  uniform  out  to  f  =  y  R/4  and  is  close  to  zero  beyond  f  =  (2-y)  R/4.  We  have  here  an 
unusual  result  as  discussed  in  Section  IV  relative  to  Figures  4-9  and  4-10.  A  far-field 
pattern  has  been  produced  in  Figure  B-l(b)  which  has  a  uniform  gain  over  an  extended  region 
in  the  mainlobe  and  a  controlled-level  stop  band  in  the  sidelobe  region.  This  unusual  charac¬ 
teristic  will  now  be  put  to  good  use. 

Let  the  array  of  Figure  B-l(a)  be  steered  to  f0  =  R/4  so  that  the  pattern  of  Figure  B-l(c) 
results.  Now  define  visible  space  by  placing  9  =  ±  tt/ 2  at  f  =  ±  (1+  y)  R/4.  This  forces  an  R 
value  of 


R  -  4/(1  ty)  (B-l) 

and  an  element  spacing  parameter  of 

T  A/R  --  (1+y)  A/4  (B-2) 

This  places  the  center  of  the  main  beam  at  0  ^90°,  As  0  decreases  the  gain  remains 

nearly  uniform  until  the  start  of  the  transition  region  at 
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is  encountered.  The  center  of  the  transition  region  is  at  broadside  (6  =  0°)  and  the  array 
gain  is  -6  dB  here.  The  stop  band  region  beings  at  9  =  and  continues  to  6  =  -90c. 

Consider  now  a  planar  array  having  line  arrays  of  the  above  type  as  elements  with  the 
axis  of  the  lines  normal  to  the  plane.  The  resulting  three-dimensional  array  does  not  require 
a  backscreen  as  the  element  pattern  properties  provide  this  function  automatically.  In  fact, 
a  polarity  change  of  the  excitation  of  the  center  element  of  Figure  B-l(a)  will  reverse  the 
element  pattern  so  that  the  main  beam  appears  at  0  =  -90°.  Thus  a  transparent  planar  array 
can  cover  two  sectors.  The  axial  symmetry  of  the  line  array  element  pattern  permits  planar 
array  steering  in  both  azimuth  and  elevation.  Because  of  reduction  in  effective  aperture  as 
the  planar  array  is  steered  toward  planar  end  fire  (which  is  at  6  =  0  in  Figure  B-l(c),  the  scan 
angle  will  probably  be  limited  by  this  factor.  This  eases  the  line  element  design  as  will  be 
seen.  As  a  practical  matter  it  would  appear  that  a  pair  of  orthogonal  vertical  planes  could 
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a)  Array  Geometry 


Figure  B-l.  End-Fire  Array  Element  Details 


easily  provide  full  azimuthal  and  reasonable  elevation  coverage.  A  third  horizontal  plane 
would  give  full  spherical  coverage. 

The  phase  taper  required  for  steering  the  line  array  to  f^  =  R/4  is  given  by  the  weight 
multipliers  w^ 

-J27rk  ~~  ~j  |  k 

w^  =  e  =  e  (B-4) 

where  the  index  k  is  now  referenced  to  the  array  end  with  w^ relating  to  the  element  at  the  left 
end  of  Figure  B-l(a).  Since  R  is  based  on  the  spacing  T,  and  since  all  of  the  even-index 
elements  are  missing  except  for  the  center  element,  all  of  the  element  currents  are  pure 
imaginary  except  for  the  center  current  which  is  entirely  real.  The  current  distribution 
shows  an  odd  symmetry  about  the  center  element.  Quite  obviously  one  can  shift  the  phase 
of  each  current  by  90°  so  that  all  currents  become  real  except  for  the  center  one  which  now 
becomes  imaginary.  Also  note  from  Equation  (B-4)  that  steering  to  f^  =  -R/4  may  be 
accomplished  by  changing  the  center-element  excitation  from  +1  to  -1.  The  bi-phase  phil¬ 
osophy  proposed  in  Section  VI  for  main-array  steering  can  apparently  provide  a  useful  func¬ 
tion  at  the  element  level  also. 

Some  examples  of  end-fire  array  designs  are  shown  in  Figures  B-2,  B-3  and  B-4.  In 
Figure  B-2  a  y  value  of  0.  3  was  used  which  required  a  T  value  of  0.  325A.  The  transition 
angle  0^  is  33°  which  would  allow  a  main-array  scan  angle  of  ±57°  from  broadside  on  each  of 
the  two  planar  faces.  This  design  would  be  suitable  for  a  main  pattern  sidelobe  specification 
of  about  -55  dB.  If  higher  sidelobe  levels  are  involved  a  simpler  end-fire  design  results. 

For  example  the  5  active  elements  involved  in  the  array  of  Figure  B-2  may  be  reduced  to  3 
for  a  -30  dB  sidelobe  level  for  the  same  scan  angle. 

Figure  B-3  shows  a  nominal  -30  dB  design  for  y  =  0.6  which  reduces  the  transition 
angle  to  15°  permitting  scanning  to  ±75°.  The  number  of  active  elements  is  still  five  which 
indicates  that  sidelobe  level  was  traded  for  increased  scan  angle  for  a  fixed  end-fire  element 
complexity. 

This  technique  can  easily  be  pushed  to  impractical  limits  as  shown  in  Figure  B-4. 

Here  y  -  0.9  which  yields  a  transition  angle  of  3°  and  a  scan  angle  of  ±87°.  This  end-fire 
design  required  39-active  elements  to  achieve  the  scan  angle  and  sidelobe  performance 
shown.  It  turns  out  that  sidelobe  level  is  not  nearly  as  expensive  as  is  scan  angle  as  regards 
active  element  count.  This  is  fortunate  since  loss  of  effective  aperture  will  limit  the  scan 
angle  in  most  practical  designs.  At  the  scan  angles  anticipated,  the  end-fire  element  designs 
appear  to  be  quite  manageable. 
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Pattern,  5  Active  Elements,  (y  =  0. 3,  T  =  0. 325A,  Scan  Angle  =  ±57°) 
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Figure  B-3.  Line  Array  Pattern,  5  Active  Elements, (y  =  0.6,  T  =  0.40A,  Scan  Angle  =  ±75°) 
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Figure  B-4.  Line  Array  Pattern,  39  Active  Elements,  (y  =  0. 9,  T  =  0.475A,  Scan  Angle  =±87°) 


